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Abstract. We present new J = 1 — ► 12 CO observa- 
j tions along the northern spiral arm of the grand-design 
spiral galaxy M100 (NGC 4321), and study the distribu- 
tion of molecular hydrogen as derived from these observa- 
tions, comparing the new data with a set of data points 
on the southern arm published previously. We compare 
these measurements on both spiral arms and on the inter- 
arm regions with observations of the atomic and ionized 
hydrogen components. 

We determine massive star formation efficiency param- 
eters, defined as the ratio of Ha luminosity to total gas 
mass, along the arms and compare the values to those in 
the interarm regions adjacent to the arms. We find that 
these parameters in the spiral arms are on average a fac- 
tor of 3 higher than outside the arms, a clear indication 
of triggering of the star formation in the spiral arms. We 
I discuss possible mechanisms for this triggering, and con- 
P elude that a density wave system is probably responsible 
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for it. We discuss several possible systematical effects in 
5j some detail, and infer that the conclusions on triggering 
are sound. We specifically discuss the possible effects of 
extinction in Ha, or a non-standard CO to H2 conversion 
factor (A), and find that our conclusions on the enhance- 
ment of the efficiencies in the arms are reinforced rather 
than weakened by these considerations. 

A simple star forming scheme involving threshold den- 
sities for gravitational collapse is discussed for NGC 4321, 
and comparison is made with M51. We find that the gas 
between the arms is generally stable against gravitational 
collapse whereas the gas in the arms is not, possibly lead- 
ing to the observed arm-interarm differences in efficiency, 
but also note that these results, unlike the others obtained 
in this paper, do depend critically on the assumed value 
for the conversion factor. 
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1. Introduction 

Although it is clear that stars are forming at higher rates 
in spiral arms than in the rest of the discs of spiral galax- 
ies, simple comparison of these rates does not tell too 
much about the physics of the massive star formation (SF) 
process. It is not clear whether the enhanced star forma- 
tion rate (SFR) in the arms is due exclusively to and de- 
pends linearly on an enhanced neutral gas density there, or 
whether specific amplifying mechanisms are at work (see 
Elmegreen 1995 for a review). One can formulate a massive 
star formation efficiency (MSFE) parameter as the SFR 
per unit gas mass; or the Ha luminosity, _L(Ha), divided 
by the surface density of the total neutral gas, molecular 
plus atomic, c(Hl+H2). Assuming the Ha is emitted in 
ionization-bounded Hll regions, _L(Ha) is monotonically 
dependent on the number of ionizing photons, and thus 
the number of ionizing stars per unit area of galactic disc, 
while c(Hl+H2) is a measure of the amount of neutral gas 
available for SF. 

SF efficiency parameters in spiral discs defined in sim- 
ilar ways, have been studied by a number of authors (e.g. 
DeGioia-Eastwood et al. 1984; Thronson et al. 1989; Wik- 
lind & Henkel 1989; Devereux & Young 1991), but for the 
problem stated above it is not satisfactory to determine 
global values for efficiencies. Tacconi & Young (1986) and 
Lord & Young (1990) discussed efficiencies in arm and 
interarm regions separately, but the resolution of their 
data was not high enough to sample the arms adequately. 
Cepa k Beckman (1990) and Knapen et al. (1992) did 
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use ionized and neutral gas data at high enough resolu- 
tion (~ 15") to cleanly map the arm regions and inter- 
arm regions separately. They defined the ratio, e, of the 
MSFE parameter in an arm, and in the neighbouring in- 
terarm disc. If the efficiencies in and outside the arms are 
the same, e equals unity (this can of course be the case 
even when the SFR in the arms is enhanced, if the un- 
derlying gas density is similarly enhanced). A value of e 
significantly higher than unity, however, implies directly 
that some mechanism amplifies, or triggers, the efficiency 
of massive SF in the arms as compared to the adjacent 
interarm region, in such a way that in the arms there is a 
non-linear dependence of the massive SFR on the under- 
lying neutral gas density. Cepa & Beckman (1990) found 
high values for e in the grand-design spirals NGC 3992 
and NGC 628, implying triggering of massive SF in the 
spiral arms. They interpreted the strongly modulated and 
symmetrical behaviour of e along the arms as signatures 
of a density- wave. 

Knapen et al. (1992) applied the technique to M51, 
including the H2 column density (which dominates within 
the inner 5 kpc) in the analysis. Results similar to those 
in the previous objects were found in M51, i.e. e showed 
peaks ranging in some cases to well over 10 in amplitude, 
whose radial distances from the centre of M51 were essen- 
tially the same for both arms, as were the distances of the 
troughs between the peaks. Here again resonance structure 
in e, as well as triggering of massive SF by the dynami- 
cal effects within the resonance pattern, were picked out 
clearly. 

Given the symmetry of the H2 molecule, which causes 
its dipolar rotational transitions to be suppressed, mea- 
surement of CO emission is the only practically applicable 
method for determining the amount of molecular hydrogen 
in the disc of a galaxy. It is measured in the emission lines 
of the CO molecule (in this case the line corresponding 
to the 12 CO J = 1 — ► transition). There are problems 
however with this method of determining H2 parameters. 
One is that the line is almost always optically thick, the 
second that the amount of observed CO may not be mono- 
tonically related to an H2 column density. We will discuss 
possible implications of these uncertainties in some more 
detail below, but will assume for the moment that we can 
convert the observed CO emission into a column density of 
molecular hydrogen using a single conversion factor (here- 
after X), which is a satisfactory approximation based on 
direct separate measurements of CO and H2 within the 
Galaxy. A review on the use of CO to measure extragalac- 
tic H2 densities is given by Young & Scoville (1991). 

In the present paper we add another galaxy to the as 
yet very limited number of objects that has been studied 
in detail with the arm/interarm MSFE technique, using 
new observational data on the ionized, atomic and molec- 
ular hydrogen for NGC 4321 (=M100). This is a grand- 
design late-type galaxy with a weak stellar bar (Pierce 
1986; Knapen et al. 1993). Earlier CO observations of this 



galaxy were made by Young & Scoville (1982), Young, 
Scoville and Brady (1985; only the central position) and 
Kenney & Young (1988) using the FCRAO antenna with 
a resolution of 50". Recently, Garcia-Burillo, Sempere & 
Combes (1994) used the 30m IRAM telescope to map a re- 
gion along the minor axis of NGC 4321 in the 12 CO 1 — >■ 
and 2 — *■ 1 transitions at 21" and 12" resolution, and used 
their observational results as input for numerical model- 
ing of the galaxy, leading to an estimate of the pattern 
speed of the density wave pattern in the galaxy (Sempere 
et al. 1995). Rand (1995) obtained interferometric CO ob- 
servations of the central region and the Southern spiral 
arm using BIMA. We present new 12 CO 1 — *■ obser- 
vations at 15" resolution with the Nobeyama millimetre 
dish of a number of points along the two main spiral arms 
of NGC 4321, and of points in the neighbouring interarm 
region (see also Cepa et al. 1992). 

In this Paper, we will denominate as the north arm the 
arm originating at the end of the bar east of the nucleus, 
passing north of the centre and continuing toward the west 
and later (but outside the range of our CO observations 
and thus of our efficiency study) the south-west. The south 
arm starts west of the nucleus at the end of the bar, passes 
the central region on the south, and continues toward the 
east and north-east. We assume a distance of 17.1 Mpc to 
NGC 4321 (Freedman et al. 1994), 1.0 kpc in the disc of 
the galaxy thus corresponds to 12". 

We will first (Section 2) present the new CO obser- 
vations, and briefly describe the H I and Ha data used. 
The CO emission is discussed in some more detail in Sec- 
tion 3, and the combination of the CO data with the Ha 
and H I data, including a determination of the efficiencies 
and their arm/interarm ratios, is described in Section 4. 
The efficiency method is discussed critically in Section 5, 
and in terms of density waves in Section 6. Section 7 dis- 
cusses a simple scheme for the observed efficiencies based 
on gravitational instabilities. The main conclusions are 
briefly summarized in Section 8. 

2. Observational Material 
2.1. CO Observations 

We observed NGC 4321 twice using the 45m telescope at 
the Nobeyama Radio Observatory (NRO) in Japan, from 
March 4 th to 8 th , 1991 and from January 18 th to 22 nd , 
1992. During the first run we observed points along the 
southern arm, as described by Cepa et al. (1992). In the 
second run we concentrated on the north arm. We present 
those new data points here. Since apart from the receiver 
used, the observational procedure was the same in both 
runs, we refer the reader to Cepa et al. (1992) for further 
details. The receiver used for the new set of observations 
was a SIS with a 2048 channel acousto-optical spectrome- 
ter providing an individual channel width of 250 kHz and 
a frequency coverage of 250 MHz. We observed the 12 CO 
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Fig. 1. Points in the disc of NGC 432f observed with the 45 m NRO telescope, plotted on a grid indicating the distance in 
right ascension and declination from the centre of the galaxy, in arcsec (centre is at RA= 12 22 m 22:8 (f 950), 8 = +f 6° 06' 00" 
(1950)). Large circles correspond to beam size (15" FWHM). Filled symbols indicate arm points, open symbols interarm points. 
The dotted line is the position angle of the major axis, PA = 153°. Numbers assigned to observed points are indicated. 
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Fig. 2. Baseline-subtracted spectra of the observed points along the northern arm of NGC 432f , binned in intervals of 5 km/s. 
Antenna temperatures are not corrected for main beam efficiency (a factor of 2.2). Point numbers are indicated in the top right 
corner of each spectrum. Interarm points are in the right-hand column and their ordinate axes are to the right. Arm points are 
in the remaining columns and their ordinate axes are to the left. 



(J= 1 — *■ 0) transition near 115 GHz. At this frequency, 
the beam size (FWHM) is 15", or slightly smaller than 
a typical arm width in this galaxy. The main beam effi- 
ciency is 0.45. Pointing was checked every ~1.5 hours us- 
ing maser and point-like continuum sources. Atmospheric 
water vapour and wind speed were low during the obser- 
vations. Final integration times ranged from 15 to 58 min- 
utes, made up from individual 20 s long integrations, with 
~50% of the points having total integrations larger than 
40 minutes. The individual spectra were co-added, and 
baseline subtracted by polynomial fitting. We determined 
total integrated antenna temperatures (J T*8v) and the 
noise therein (r.m.s. of line-free channels). The observed 
points, with our assigned reference numbers, are shown 
in Fig. 1, in an RA-dec grid centred on the nucleus of 
the galaxy. The baseline subtracted spectra are shown in 
Fig. 2. 



In the rest of this paper (unless specifically indicated), 
we use a constant CO-H2 conversion factor X of 2.8 x 10 20 



[K km 



_1 (Bloemen et al. 1986). The calibration 



(from / T*6v to N(H j) is then 

N(H) = 2.8 • 10 20 x 2 x — x / T*6v, 
Vi J 

with N(H) in atomscm -2 . The factor 2 in the calibration 
is required because the H2 molecule consists of two hydro- 
gen atoms, and rji, is the main beam efficiency, rji, = 0.45. 

The results are summarized in Table 1 (note that 
points 0-18 were published in Cepa et al. 1992). The first 
column is the number assigned to a given sky point. The 
distance to the centre r (in arcsec) and the position an- 
gle 6 on the plane of the sky (measured from N towards 
E) are given in columns 2 and 3. Column 4 is the depro- 
jected distance to the centre ro in seconds of arc. We used 
an inclination angle i = 27° (de Vaucouleurs et al. 1976) 



J. H. Knapen et al.: Molecular gas observations and enhanced massive star formation efficiencies in MfOO 5 



and a major axis position angle of 153° (Knapen et al. 
1993) for the deprojection. Columns 5 and 6 are the offsets 
from the centre (at position RA= 12 h 22 m 22!8 (1950), 6 = 
+ 16° 06' 00" (1950)) in right ascension and declination, in 
seconds of arc (north and east are positive). Column 7 is 
the integrated antenna temperature J T*8v, and its error 
(the errors were evaluated from the line-free channels). 
T* is not corrected for the main beam efficiency. The last 
column (8) shows whether the points are measured on an 
arm or in the interarm disc; A means arm point, I interarm 
point, and C the centre of the galaxy. 

2.2. Combination with Ha and Hi data 

The H I data used here were described in detail by Knapen 
et al. (1993). We use here the VLA total Hi map at a 
resolution of 15", comparable to the resolution of the CO 
data. The map units of mJy/Bxm/s can be transformed 
into N(H) by using the conversion 

1 mJy/B x m/s = 4.89 x 10 18 atoms cm" 2 . 

We used a new Ha continuum subtracted image of 
NGC 4321 obtained with the TAURUS camera in imag- 
ing mode on the 4.2m William Herschel Telescope. The 
observations and data reduction are described in detail 
by Knapen (1992) and Knapen (1996). In Fig. 3, we show 
a contour representation of the Ha map of NGC 4321. 
One instrumental count in the Ha image corresponds to 
2.46xl0 33 erg/s. 

We now describe the additional steps followed to com- 
bine the CO, H I and Ha data. In the case of CO, we have 
only a limited number of data points, whereas in the case 
of Hi and Ha, we can use a map covering the whole disc 
of NGC 4321. For a direct comparison of the three data 
sets, we are thus limited to the CO positions. The common 
spatial resolution used for all observations is ~ 15". Since 
the Ha map has a considerably higher intrinsic resolution 
(~ 1"), we convolved the whole map to obtain a 15" reso- 
lution image, using a gaussian "seeing" profile. In so doing 
we have degraded the resolution, and thereby the intensity 
of the Ha peaks and the arm/interarm contrast, to make 
the Ha directly comparable to the H I and CO data. Note 
that combining data sets of different resolutions will yield 
unphysical results. 

The Hi and Ha images, now both at 15" resolution, 
were transformed to a common pixel grid. We chose to 
use 512 x 512 pixels of 1" x 1" for the combined grid, 
where the pixel sizes of the original maps were 4" x 4" in 
the case of the Hi map, and 0"27 x 0"27 in the case of 
the Ha image. Using positions of field stars, we first gave 
the maps a common orientation and centre-point. Subse- 
quently, the re-gridding was done using a bicubic interpo- 
lation algorithm. The Ha luminosity and the H I column 
density were determined directly from the re-gridded H I 
and Ha images. Average values over small boxes centred 
at the positions observed in CO were taken as the H I and 



Ha on-arm points. The rms of the > 100 pixels over which 
they were determined was taken as an indication of the er- 
ror value in the points. A value for the interarm Ha or H I 
intensity for a specific point was determined as an average 
of measured intensities in the interarm regions measured 
at both sides of the arm. In most cases, values determined 
at opposite sides of the arms did not differ much. In the 
case of the second half of the north arm (the part toward 
the NW, or at position angle along the arm of > 60°), 
the determination of reliable interarm values was difficult, 
especially in H I. This is caused by the fact that the arm is 
not well defined there. In the case of Ha, it was still pos- 
sible to determine the interarm emission for each point 
separately. In the case of Hi though, we assumed a fixed 
interarm value for the complete second part of the arm, 
determined as an average value over that whole interarm 
region. 

3. CO Emission 

3.1. Radial behaviour 
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Fig. 4. a. (top) CO integrated intensity (J T*8v, in Kkm/s) 
as a function of deprojected distance to the centre m of 
NGC 4321. Arm points are filled triangles, interarm points 
open triangles, b. (bottom) As Fig. 4a, now for log ( j T*8v) 
as a function of rr>- Dotted, dashed and solid lines are expo- 
nential fits, with h = 46", h = 59" and h = 76", respectively 
(see text). 

In Fig. 4a we show the observed CO intensity as a function 
of deprojected radial distance for all measured points. In 
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Fig. 3. Contour plot of the part of the continuum subtracted Ha image of NGC 432f used in the present paper, shown at the 
same scale and orientation as Fig. f. Contour levels are 0.5, f.5, 4.4, f3.3 and 26.6 x 10 36 ergs -1 . The straight feature near 
RA= +80", 8 = —60" is an artifact in the image (probably a meteorite track). 



Fig. 4b we plot the same points as in Fig. 4a, but now 
on an exponential scale, indicating exponential fits with 
straight lines here: the dotted line represents a fit with the 
central point included (hco = 44 ± 8"), the dashed line a 
fit where the centre was excluded (hco = 57 ± 12"), and 
the solid line a fit to all points excluding the central point 
and those in the bar (hco = 76 ± 14"). The values for the 
first two cases are only marginally different from the values 
given by Cepa et al. (1992) and Kenney & Young (1988). 
Young & Scoville (1982) determined a scale length of 58" 
fitting their 3 disc points. Our value is also consistent with 
the minor axis CO profile shown by Garcia-Burillo et al. 
(1994), although they do not explicitly fit an exponential 
to their data. Nakai (1992) described radial CO profiles 
in some other barred spirals, showing that the behaviour 
within and outside the bar region is different. Our fit to 
the points outside the bar indicates that indeed the profile 
is steeper when including the bar points. The presence of 
a distinct central component is clear, but we claim this 
is a (weak) bar component, and not just a central peak 
as claimed by Kenney & Young (1988) and Cepa et al. 
(1992). 

3.2. CO intensity along the arms 

In Fig. 5 we plot the integrated CO intensity against the 
position angle along the arms. The two main arms, and 
arm/interarm points are shown in different point styles. 
The arm points are connected by continuous lines, the 
interarm points by dashed lines. A position angle of 0° 
in the figure corresponds to the starting point of the two 
main spiral arms, defined here at the end of the bar (our 
CO data points 3 and 23). This zero-point corresponds to 
an angle of 120° (for the north arm), and 300° (south arm) 
on the galaxy (measured north over east). The position 
angle along the arm is measured clockwise. 

It is immediately clear that the behaviour shown in 
Fig. 5 is quite different from the exponential fall-off seen 




o 
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Position angle along arms 

Fig. 5. CO integrated intensity (J Ta8v, in Kkm/s) along the 
two main arms of NGC 4321. Filled symbols are arm points, 
open symbols are interarm points along the same arm. The 
points along the south arm are represented by circles, points 
along the north arm by triangles. 



in the previous section. The CO emission shows variations 
along the arms, but no generally falling trend is seen. This 
is mainly due to the fact that our arm points are con- 
centrated in a small range of galactocentric radius. Al- 
though the sampling is limited, the interarm CO emission 
looks rather constant along the two arms. To estimate 
the interarm emission at a given point along an arm, de- 
termined here as the average between interarm emission 
from both sides of the arm point in question, it is im- 
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portant to keep radial distance effects in mind. Along the 
north arm, we measured two pairs of two interarm points 
on each side of the spiral arm. In Fig. 5 it is seen that 
f T*8v = 2.0 ± O.lKkm/s for both the points at larger 
distances than the arm, while J T*8v = 4.1 ± 0.3Kkm/s 
for the points closer to the centre. We will use an average 
interarm CO emission of J T*8v = 3 ± 0.5 Kkm/s for the 
whole north arm. Although the sampling is limited, this 
is justified by the symmetric distribution of the interarm 
points and by the fact that the emission from the first set 
of interarm points is almost equal to that from the second 
set. In the case of the south arm, all the interarm points 
measured lie somewhat further out in the disc than the 
arm points. We assume a general interarm value for the 
south arm of f T*8v = 3 ± 0.5 Kkm/s, the same as that 
found for the north arm. 

The profiles along the two arms show some pronounced 
differences. Along the south arm, the CO emission shows 
two broad peaks, centred at position angles of 30° and 
120°, and reaching peak values of some 2 — 3 times the local 
interarm emission (Fig. 5). At either side of these peaks, 
the emission falls off to a level of f T*8v = 3 — 4 Kkm/s, 
equal to or only marginally higher than the interarm emis- 
sion near those positions. The CO emission along the 
north arm looks quite different from that along the south 
arm. It is almost completely constant at a level of some 
3 — 4 Kkm/s, but shows two narrow peaks, near position 
angles 0° and 60° . Note however that both of these peaks 
consist of only one observed point. 

The arm/interarm ratio varies between 1 and 3. 
This value is lower, but only a little lower, than the 
arm/interarm ratio in CO of ~ 3 as measured in M51 (see 
Knapen et al. 1992). We will discuss the arm and interarm 
CO profiles in more detail in the next Section, where we 
will combine the CO data with Hi and Ha observations. 

4. Gas distribution and efficiencies along the arms 

4.1. Ha, Hi and H^ along the arms 

The arm and interarm gas densities and Ha luminosities 
at the positions of the points observed in CO are shown 
in Fig. 6a for the south arm, and in Fig. 6b for the north 
arm, plotted as functions of position angle along the arms. 

We will first discuss the south arm. The H2 and the 
Ha show a similar structure along this arm. A depression 
is seen between some 70° and 100° in position angle; fur- 
ther along the arm (higher PA) both Ln a and <7h 2 are 
enhanced. In the first part of the arm, Ln a is slightly 
enhanced, and the <7h 2 possibly as well. The atomic gas 
density is in general a factor 10 lower than the molecular 
gas density (this value depends on the exact value of X). 
Thus, as in the case of M51, the role of variations of uhi in 
the MSFE will be limited. The uhi along the arm shows 
hints of similar behaviour to the <7h 2 and Ln a , in partic- 
ular it shows a depression near PA= 80° which coincides 



with the depression seen in the ionized and molecular gas 
tracers. There is also a depression in uhi near PA= 20°, 
where we noticed the low CO emission discussed before. 
The interarm values vary little with increasing position 
angle along the arms. Arm/interarm contrasts are gener- 
ally about 7 in Ha, about 2.5 in Hi and about 2 in H2. 

The distribution of the gas tracers in the north arm 
(shown in Fig. 6b) is quite different from that in the south 
arm. The Ha luminosity is seen to start at relatively high 
values when beginning at PA=0° along the arm. This is 
a region at the end of the bar where the SF is enhanced. 
Outside this region, Ln a falls off to low values, and stays 
low until about PA=150°. In this whole part, the north 
arm is poorly defined: the Ha arm/interarm contrast is 
only between 1 and 3. At the end of the region observed in 
CO, the Ha luminosity is larger. This is the beginning of a 
part of the arm west of the nucleus, which is strongly emit- 
ting in Ha. The H I density shows little structure along the 
north arm, and is generally at a lower level than the south 
arm. The Hi arm/interarm contrast is only about 1.5 in 
the first part of the arm, and drops even further, to values 
close to unity, in the middle of the arm region observed in 
CO and considered here, uhi seems to rise together with 
Ln a towards the end of the arm. As discussed before, the 
interarm emission is hard to determine in the outer part of 
the arm region considered in the present study. For those 
outer points we have assumed an average value over that 
whole part of the interarm region as a generic interarm 
value for <thi- The molecular hydrogen density is gener- 
ally hardly enhanced with respect to the interarm value 
along the north arm. In general, we can remark that the 
north arm of NGC 4321 does not seem very interesting 
in terms of enhanced SF. It is difficult to define its form 
exactly due to the low contrast with the surrounding in- 
terarm regions. As can be seen from Hi and particularly 
Ha maps, the massive SF activity begins to rise where our 
CO observations had to stop. 

4-2. MSFE arm /interarm ratios along the arms 

Following Cepa & Beckman (1990) and Knapen et al. 
(1992), we define a massive star formation efficiency 
(MSFE) parameter as the Ha luminosity divided by the 
total (atomic plus molecular) gas density at a certain point 
in a galaxy, which is closely proportional to the massive 
SFR per unit underlying gas mass. The arm/interarm ra- 
tio of the MSFE, e, is then defined as 

_ MSFE^ _ (Lu a ) A (a hi + <th 2 )ia ,^ 
MSFV IA ~ (L Ha ) IA X (a HI + (Th 2 ) A ' 1 ' 

where the subscripts A and I A stand for arm and in- 
terarm, respectively. The values for e along the two arms 
are plotted as a function of position angle along the two 
arms of NGC 4321 in Fig. 7. 

The dashed line in the figure, at an e of unity, indicates 
the case of equal MSFE's inside and outside an arm. This 
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Fig. 6. a. Ha (top panel), Hi (middle panel) and H2 (lower panel) distributions along the south arm of NGC 4321. For Ha 
and H I, filled symbols are on-arm values, and open symbols indicate corresponding interarm values. The interarm value in H2 
is indicated by the dotted line in the lower panel. Note that values along the ordinate are an order of magnitude higher in H2 
than in Hi. b. As Fig. 6a, now along the north arm of NGC 4321. 



implies that if the massive SFR is higher in the arms, the 
underlying gas density is proportionally higher. As can be 
seen in Fig. 7, the values off, as determined from the Ha, 
H I and H2 data points described in the previous para- 
graph, are almost always significantly higher than unity, 
and in no case significantly lower. This indicates that there 
is more massive SF in the arms than might be expected 
from the amount of neutral gas available, as compared to 
the interarm regions. It thus means that the massive SF 
in the arms is enhanced, which implies that some specific 
triggering process is at work there. 

An average value for e in the south arm is e = 3.5±0.5, 
whereas in the north arm it is somewhat lower, e = 
2.5 ± 0.5. It is important to discuss in some detail the 
peaks and dips seen in e along the arms, since such peaks 
could originate in less reliable determinations of one of the 
gas tracers. An example of such a case may be the peak 
seen in the south arm at some 20° in position angle. This 
peak is not a result of enhanced massive SF in the arm, 
but rather of the low CO emission observed at that point 
(detection with low S/N). In the north arm, two sepa- 



rate points that were observed in CO showed significantly 
higher emission than their neighbour points. The result of 
the one-point CO peaks in the e run along the north arm 
is a pair of relative dips, at positions angles of some 10° 
and 70° along the arm. Finally, there is the peak in e seen 
at the eastern end of the north arm (at PA~ 170°). At this 
position, Ln a is higher than at previous points, because 
one is starting to see the bright star forming complex west 
of the nucleus. Efficiencies may well be enhanced in this 
region, but more CO observations further along the arm 
are needed to confirm this. 

Keeping these uncertainties in mind, the relations 
shown in Fig. 7 can be summarized as follows. In both 
arms, the values of e are clearly and significantly higher 
than unity. Typical values for the arm/interarm contrast 
in MSFE in the part of the disc of NGC 4321 studied 
here are around 3. There is no convincing evidence for 
the existence of peaks and dips in e along the arms of the 
high amplitudes seen in M51 (of order 10-20; Knapen et 
al. 1992), or even of lower amplitude. No symmetric pat- 
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Fig. 7. Arm/interarm ratio of the MSFE (e) along the two 
main arms of NGC 4321: north arm in the upper panel, south 
arm in the lower panel. Dashed line at e = 1 indicates equal 
efficiencies in and outside the arms. 

terns, and in fact no clear patterns at all, are seen in the 
runs of e for the two different arms. 

Equation (1) shows that high values of the 
arm/interarm ratio of the MSFE can in principle be 
caused by any of the four factors: Ha luminosity and to- 
tal gas density, in and outside the arms. Enhanced Ln a on 
the arms or <7 gas between the arms, or reduced interarm 
Ln a or on-arm <7 gas will produce higher e's. A combina- 
tion of more than one of these factors can reduce or am- 
plify the effect. For NGC 4321, the profiles in Fig. 6 show 
that the enhanced efficiency in the arms is due only to an 
arm/interarm contrast in Ha higher than in neutral gas 
(thus again to enhanced on-arm Ha emission). Interarm 
Ha emission or neutral gas densities do not "conspire" 
to suppress any Ha peaks which might in principle be 
present, but which do not show up in the run of e shown 
in Fig. 7. The values of e larger than unity do indicate an 
enhanced MSFE in the spiral arms. 

5. Discussion: on the validity of e as an indicator 
of enhanced on-arm efficiencies 

5.1. The CO to H2 conversion factor X 

The use of arm/interarm ratios in e mitigates problems 
due to the need to convert observed CO column densities 
to H2 column densities. While there is evidence that X 
does not in fact vary by large factors across most of the 
disc of the Galaxy (Bloemen et al. 1986; Scoville & Good 



1989) , and that the Galactic value is valid in some external 
galaxies (Vogel, Boulanger & Ball 1987; Wilson & Scoville 

1990) , there are also indications that X may vary within 
the Galaxy (e.g. Polk et al. 1988; Digel et al. 1995) or 
among galaxies (e.g. by as much as a factor of three in the 
well-studied galaxy M51: Adler et al. 1992; Rand 1993a; 
Nakai & Kuno 1995). There is however no evidence to date 
that X is different in NGC 4321 (Rand 1995). There is no 
evidence for, nor against, arm-interarm variations in X 
in general. The fact that e as used in the present paper 
contains a ratio of <7H 2 (arm)/<7H 2 (interarm), with both a 
at the same galactocentric radius, certainly means that a 
different value for X for NGC 4321, or a radial gradient 
in X , will have no effect on our determination of e. 

A potentially more fundamental problem with the use 
of CO as an indicator of the amount of molecular hydrogen 
has been exposed recently in a number of papers. From a 
study of molecular clouds in the SMC, Rubio, Lequeux 
& Boulanger (1993) concluded that X is significantly dif- 
ferent there from the "standard" Galactic value, and that 
the CO luminosity depends primarily on the UV radiation 
field and, to a lesser extent, on the cosmic-ray (CR) flux 
(see also Lequeux et al. 1994). Allen & Lequeux (1993) de- 
tected CO in the inner disk of M31. The UV radiation field 
and the CR flux are both very low in the regions where 
the "cold" CO was observed, and kinematic temperatures 
of the clouds can drop to less than 5 K (Allen et al. 1995). 
All this obviously implies that X changes with the UV 
and CR flux, with in general relatively more CO luminos- 
ity per unit molecular hydrogen in zones with higher UV 
and CR flux. Note that X may increase again in regions of 
very high UV flux, due to the higher degree of destruction 
of CO compared to that of H2, and that also abundances 
may affect the CO luminosity. 

Although the extragalactic observational evidence for 
these ideas is as yet limited to the SMC and M31, it is 
worthwile to discuss their possible influence on the results 
of the present paper. Our main conclusion on triggering of 
the massive SF in the arm regions is based on our obser- 
vations of an arm/interarm contrast that is larger in Ha 
than in CO. We can expect that both the UV radiation 
field and the cosmic-ray flux are larger in the arms than 
between the arms, and thus the CO luminosity would be 
higher in the arms, following the reasoning laid out in the 
previous paragraph. In that case, X would be higher in 
the interarm regions, and there would be more molecular 
hydrogen between the arms than one would assume using 
a constant X. But this means that the arm/interarm con- 
trast in molecular hydrogen is lower than assumed, that 
e is higher, and thus that the values for the arm/interarm 
efficiency as determined in Section 4 are in fact lower lim- 
its to the true values. This reinforces our conclusions on 
triggering of the massive SF in the arms. 

The fact that in the outer parts of the disc of 
NGC 4321 (and also M51), where the neutral gas den- 
sity is dropping below 10 21 cm -2 , and where Hi begins 
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to compete with H2 for dominance in column density, the 
behaviour of e shows values consistently in the range close 
to 3, indicates that this value for e cannot be an artifact of 
a wrongly adopted value for X . We feel it safe to conclude 
that our evidence for the triggering of the massive SF in 
the arms is sound. 

5.2. Extinction by dust 

Dust extinction could in principle affect our measurements 
of the Ha emission in the arm and interarm regions, in 
the sense that one would expect more dust in the arms, 
and thus underestimate the arm Ha emission compared to 
the interarm emission. If this were the case, the intrinsic 
arm/interarm ratio in Ha would be larger than measured, 
and this means that the values for e can only be larger than 
determined here. Beckman et al. (1995) studied in detail 
the dust content and distribution in NGC 4321 through 
the comparison of scale lengths in the optical, separately 
for arm and interarm regions. They found that there is 
not much dust in this galaxy, that there is significantly 
less dust between the arms than in them, and that the 
dust in the arms is mostly confined to narrow dust lanes. 
It thus seems safe to assume that the interarm Ha is very 
little, if at all, affected by dust extinction. The Ha emis- 
sion in the arms arises mostly from regions outside the 
well-defined dust lanes, so the effect of extinction is also 
small in the arms. However, if we are in fact underesti- 
mating the on-arm Ha, e.g. because of emission arising 
inside the dust lanes which will thus be heavily obscured, 
we are also underestimating the arm-interarm contrast in 
Ha, and thus the efficiency ratios. The absence of peaks, 
as seen e.g. in M51 (Knapen et al. 1992), in Ha along the 
arms of NGC 4321 cannot be caused by dust extinction, 
since this would imply the localized presence of huge quan- 
tities as dust which would show up in optical images, in 
the work of Beckman et al. (1995), and/or as strong neu- 
tral gas concentrations. There is no observational evidence 
for any of these effects, but in any case their presence 
would mean that our numbes for arm/interarm contrast 
in MSFE were lower limits and our qualitative conclusions 
would be strengthened, not weakened. 

6. Efficiencies and density waves in NGC 4321 

Comparing the results for the arm/interarm ratio of 
MSFE along the arms of NGC 4321 with those along the 
arms of M51 (Knapen et al. 1992), some important differ- 
ences are obvious, but some common features must also 
be mentioned. Whereas in M51 we observed a strong sym- 
metric pattern of high, narrow peaks and low dips (with 
values close to unity) in the run of e along the arms, in 
NGC 4321 no pattern can be seen, and e is practically 
constant along the arms. High peak values in e of 10-20 
as seen in M51 are not convincingly present in NGC 4321, 
where e takes moderate values, with the exception of one 



individual point on each arm, where the CO intensities are 
lowest, and consequently the uncertainties highest. Also, 
e does not drop down to unity in the case of NGC 4321 
(except for a few points in the north arm), whereas it does 
in M51, between the peaks of high efficiency. On the other 
hand, an average of e over the arms in M51 is similar to 
that determined for NGC 4321, of about 3 (±0.3). 

Thus in NGC 4321, as in M51, the MSFE is enhanced 
in the arms with respect to the interarm regions. We can 
conclude that some kind of triggering mechanism enhances 
the massive SF in the arms of NGC 4321 (otherwise e 
would be ~ 1), but that the resonance pattern in e is not 
apparent. Basically, what we have measured in NGC 4321 
is that the arm/interarm ratio in Ln a is larger than that 
in <th 2 , without any significant changes along the parts of 
the arms studied (note that, even though dependent on 
the exact X value, the contribution of atomic hydrogen is 
of much reduced importance here, since it is far less abun- 
dant than molecular gas). In M51 (Knapen et al. 1992), 
and also in NGC 628 and NGC 3992 (Cepa k Beckman 
1990), the enhanced efficiencies in the arms could be di- 
rectly related to the presence of a density wave organizing 
the massive SF, through the symmetric pattern of peaks 
and dips in the run of e along two main arms. The dips 
could in all these galaxies be identified with density wave 
resonances, and the global nature of this mechanism nat- 
urally explains the remarkable observed symmetry. 

In NGC 4321, where no pattern in e is seen in ei- 
ther arm, the efficiencies in the arms may still well be 
caused by density wave triggering. There is evidence from 
other observations that a density wave system is present 
in the disc of this galaxy, but also that its main resonances 
may be outside the region studied in the present paper. 
Streaming motions over the arms have been observed in 
Hi (Knapen et al. 1993). Elmegreen, Elmegreen & Seiden 
(1989) were able to fit density wave resonances to radial in- 
tensity variations in arm profiles, as determined from com- 
puter enhanced optical images, and placed the co-rotation 
radius (Rcr) at ~ 110" (well beyond the ends of the bar). 
This value for Rcr has been confirmed by Sempere et al. 
(1995), who applied a two-dimensional test for reversal of 
streaming motions (following Canzian 1993) to an H I ve- 
locity field. However, Rcr was estimated at around 75" 
by Arsenault et al. (1988), using an ionized gas rotation 
curve, and at 75" ± 7" by Canzian & Allen (1995), who 
optimized the fit of spiral density wave response functions 
from linear theory to an Ha velocity field. This indicates 
a degree of confusion about where the co-rotation radius 
is located in this galaxy. It may be outside the bar region, 
and possibly outside the region studied in the present pa- 
per. Since a double inner Lindblad resonance is located 
well inside the region studied here (Arsenault et al. 1988; 
Knapen et al. 1995), we may be plotting e between two 
resonances, and not covering the ILR or CR, which could 
explain why we do not see a clear pattern in the efficien- 
cies. 
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7. Arm vs. Interarm Gravitational Instability 

Galaxies with two symmetric spiral arms, such as 
NGC 4321 (where the general arm shape is symmetric, 
even though the distribution of SF regions in them is not) 
or M51, are believed to be controlled by spiral wavemodes 
(Bertin 1991), possibly organized by bars or companions. 
In the spiral arms, gas collapses into giant molecular com- 
plexes along an arm and may form OB stars (Scoville, 
Sanders & Clemens 1986). We see such a sweeping of gas in 
NGC 4321 as an enhancement of the neutral (molecular) 
gas density along the arms, as compared to the interarm 
regions. These clouds may lose enough of their support- 
ive energy (turbulent, rotational or magnetic) or become 
massive enough to collapse under their own gravity, and 
may then form stars. Such triggering by gravitational in- 
stabilities is described in more detail by Elmegreen (1979, 
1994) and Cowie (1981). Alternative mechanisms to in- 
crease the SF efficiency in the arms have been described 
(e.g. cloud-cloud collisions, Scoville et al. 1986). A review 
can be found in Elmegreen (1992). 

Kennicutt (1989) studied massive SF in discs of spiral 
galaxies using information on the total gas content and 
on the massive SFR. He found that the massive SFR in 
the inner, dense regions can be described by a Schmidt 
(1959) power law R = ap n g , with n = 1.3 ± 0.3, where R, 
the massive SFR, depends only on the local gas density 
p g . This relation breaks down in the outer parts of the 
galaxies Kennicutt studied, where the gas density drops 
below a critical threshold value. Massive SF is suppressed 
at densities well below this threshold, as the gas is stable 
against gravitational collapse at these densities. 

An attractively simple explanation for the observed 
triggering in the arms of NGC 4321, and also in those of 
M51 (Knapen et al. 1992) might be that the gas density is 
higher than the local critical value in the arms, and falls 
below the threshold in the interarm regions. We can in 
fact check this hypothesis with the data at hand. 

Critical densities <r c for several points along the arms 
and in the interarm regions can be derived following the 
formulae given by Kennicutt (1989). To make an estimate 
of a c , we assume that the rotation velocity is constant at 
v c = 220 kms -1 (from the Hi rotation curve of Knapen 
et al. 1993). This is a reasonable approximation for all the 
points we consider, which are in a small range of radius 
around 70". We then find that a c = 3.2 • 10 3 /x M pc" 2 , 
where x is the distance from the centre, in arcsec. We 
determine total gas densities from the CO and H I obser- 
vations (see above), add the densities in units of N(H), and 
multiply by 1.36 to account for the contribution of Helium 
to the total gas density, following Kennicutt (1989). 

In Table 2, we show the gas densities and critical 
(threshold) densities for several arm and interarm points 
in NGC 4321. We express the results in solar masses per 
square parsec (10 21 N(H) cm -2 = 8.00 JVI pc -2 ). In the 
south arm, the arm points have total gas densities about 



twice the critical density. For the interarm points a is just 
below the critical density in point 10, and slightly above 
it in point 18. These findings are in agreement with those 
of Rand (1995), who determined first the arm/interarm 
contrast in CO from streaming motions as observed with 
BIMA, and then surface and critical densities at arm and 
interarm points. In the north arm, which is, as we saw be- 
fore, weaker than the south arm both in massive SF and 
gas density, the arm points have a higher than a c for all 
the points but one, but only slightly higher. The interarm 
points though have gas densities that are below the criti- 
cal densities by about a factor of two. The last point along 
the north arm, point no. 39, shows a greater contrast be- 
tween a and a c . We saw before that at this point also e 
is enhanced, thus the enhanced massive SF here may be 
related to the lower critical density. 

We thus find that there is a significant difference be- 
tween arm and interarm points, where gas densities of the 
arm points are either well above (south arm) or similar 
(less developed north arm) to the critical density. Interarm 
gas densities are very close (south) or well below (north 
arm) a c . In the arms the gas would thus not be stable 
against gravitational collapse, and in the interarm regions 
it might just be, which could well be the reason for the en- 
hanced MSFE seen in the arms. This effect is more promi- 
nent along the south arm, which as we saw before shows 
substantially more SF than the north arm. 

For M51, Rand (1993b) estimated the value of the 
epicyclic frequency k in the arm and in the interarm re- 
gion, taking account of the streaming motions caused by 
the density wave. Rand (1993b) also used his results for 
k to estimate an arm/interarm contrast in CO, and com- 
bined those numbers with Lord's (1987) data to obtain a 
gas density. He found that whereas the arm gas is gener- 
ally unstable against collapse (<th 2 > c c ), the interarm gas 
density is near the critical value. Oey & Kennicutt (1990) 
show the radial dependence of critical and observed gas 
densities for M51, using Lord's (1987) data for the CO. 

As an extra test, we determine gas and critical den- 
sities for M51 by directly measuring the total gas densi- 
ties from the Hi (from Rots et al. 1990) and CO data 
(from Nakai et al. 1991) used by Knapen et al. (1992) at a 
few arm and interarm positions, and compare our results 
with those found by Oey & Kennicutt (1990) and Rand 
(1993b). The results for the three arm points (called Al, 
A2 and A3) and three interarm points (II, 12 and 13) con- 
sidered here are shown in Table 3. Points A3 and 13 are 
at a distance of 2.7 kpc in order to compare with Rand 
(1993b). Table 3 shows that the observed gas densities in 
both arm and interarm points lie well above the critical 
density. The fact that some of our points also lie well above 
the points determined by Oey & Kennicutt from Lord's 
data is probably an artifact of the lower resolution (45") 
of Lord's observations, or of differences in the calibration 
of the data. Rand (1993b) listed critical and inferred sur- 
face densities at a distance of 2.7 kpc. For the arms he 
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finds values for the gas density of 140 and 170 M pc -2 
(not included in his numbers is the contribution of atomic 
gas, which is ~10%), which are comparable with our 196 
M pc -2 . Rand gives a value of cr c =75-80 M pc -2 for 
the critical density, far lower than the measured density. 
For the interarm region, Rand (1993b) finds a = 45 M 
pc -2 and <7 C = 40 M pc -2 . Our value of a = 78 M 
pc -2 for the interarm point at R = 2.7 kpc (which is a 
typical point for the interarm region at that galactocen- 
tric distance) is clearly higher. The difference may again 
be due to the lower resolution or different calibration of 
Lord's (1987) observations that Rand (1993b) used, or to 
the method Rand used to derive the arm/interarm con- 
trast. 

Explicitly considering the impact of changing X, we 
useavalueof#(H 2 )//co = 1.1 xlO 20 cm" 2 [Kkms" 1 ]" 1 , 
determined by Nakai & Kuno (1995) from a comparison 
of visual extinctions and CO data toward 30 H II regions 
in M51, a result which is in general agreement with con- 
clusions from studies of virial and CO masses by Adler 
et al. (1992) and Rand (1993a). The observed and critical 
densities are much closer using the low value of X . Two 
of the arm points (Al and A2) continue to have a 2 — 3 
times larger than <r c , with A3 having a gas density only 
marginally larger than the local <r c , while gas densities at 
the interarm points are all very close to the local criti- 
cal densities, with a somewhat larger than a c for II, and 
somewhat smaller for 13. This is reminiscent of what we 
found for NGC 4321, where the arm gas may be unstable, 
but the interarm gas stable against gravitational collapse. 

It is clear that any conclusion which is directly depen- 
dent on absolute measurements of quantities of molecular 
gas in external galaxies, such as those described in this sec- 
tion, but also e.g. the study of Kennicutt (1989), depend 
critically on the adoption of a correct for X. As we saw 
before (Sect. 5.1) there are indications that the standard 
conversion factor may not be correct in general. In M51, 
there is now mounting evidence that X is in fact lower, 
as taken into account here. We note again, however, that 
so far there are no indications that X is different from 
the "standard" value in NGC 4321 (see Rand 1995), but 
also that our conclusions on triggering (Sect. 4) seem ro- 
bust against variations in X , unlike the results on gas vs. 
critical densities. 

8. Conclusions 

We have measured the 12 CO J= 1 — *■ emission from 
a sample of points along the two main spiral arms of 
NGC 4321, and from a number of interarm points, us- 
ing the Nobeyama 45 m millimetre wave telescope. In this 
paper we complement the data points presented by Cepa 
et al. (1992) for the south arm with new measurements 
along the north arm. The exponential decline of the ra- 
dial distribution of the CO integrated antenna temper- 
ature, reported before, is confirmed, but the inner part 



of the radial profile is determined by emission from the 
weak bar of the galaxy. Interarm emission is constant at 
a level of f T*8v = 3 ± 0.5 Kkm/s along the parts of the 
arms studied here. The south arm shows generally higher 
emission than the north arm. 

Combining the CO data with new H I and Ha maps, we 
have determined massive star formation efficiencies (de- 
fined as Ha luminosity per unit neutral gas mass) along 
the two main spiral arms, and especially the arm/interarm 
ratio e of these efficiencies. This is only the second galaxy, 
after M51 (Knapen et al. 1992) where e has been deter- 
mined along the spiral arms, including both the molecular 
and atomic hydrogen gas components in the analysis. As 
in M51, we find in NGC 4321 that e is larger than unity, 
with an average value of 3 ± 0.5, very similar to the aver- 
age value for M51. This implies that some mechanism trig- 
gers the massive star formation in the arms. Apart from a 
few localized peaks, values of e are almost constant along 
both arms, and somewhat higher in the south than in the 
north arm. A density wave system is a good candidate 
for the triggering agent, even though direct evidence such 
as a symmetric pattern in e similar to that found in M51 
(Knapen et al. 1992) is absent in NGC 4321. 

We find that in NGC 4321, measured gas densities are 
equal or larger than local critical densities for gravita- 
tional collapse in the arms, whereas the interarm gas den- 
sities are equal or lower than a c . We confirm earlier find- 
ings that in M51 both arm and interarm gas densities are 
higher than local critical densities, assuming the standard 
X-value. But assuming the lower value for X reported in 
the literature for M51, a picture emerges that is very much 
like the one we find in NGC 4321. Thus in both galaxies 
gravitational instability could be a mechanism involved in 
causing the excess of massive SF in the arms, as compared 
to the interarm regions. 

Our conclusion on triggering of the massive SF in 
the arms seems sound. It basically reflects a higher 
arm/interarm contrast in Ha than in molecular gas, which 
is dominant by an order of magnitude over atomic gas in 
the inner region of the galaxy studied here. From the pro- 
files along the arms for the several gas tracers, we see no 
evidence that the enhanced on-arm efficiencies are due to 
any effect other than enhanced on-arm massive SF. Dust 
extinction or changes in the CO to H 2 conversion factor 
(X) are among the effects we have eliminated as possible 
causes for enhanced efficiency in the arms. Both effects, if 
at all present, would increase the arm/interarm efficiency 
ratios and thus reinforce rather than limit our conclusions 
on triggering. 
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Table 1. Summary of CO Observations. Columns are (respec- 
tively) point number, galactocentric distance, position angle 
and deprojected distance from the centre of the galaxy, RA 
and dec offsets from the centre of the galaxy, CO integrated 
flux and error, and arm (A), interarm (I) or centre (C) point 
identification. 



Table 2. Surface and critical gas densities in NGC 4321. Col- 
umn 1: point identification as in Table 1 and Fig. 1. Column 2: 
point is in south (S) or north (N) arm, and (column 3) point 
is in the arm (A) or in the interarm (I) region. Column 4: de- 
projected distance from the centre, in arcsec. Column 5: total 
gas density, in Mq pc -2 . Column 6: critical or threshold gas 
density, in Mq pc -2 . 
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Table 3. Critical densities a c as determined by Oey & Ken- 
nicutt (1990) and Rand (1993b; for points at 2.7 kpc from the 
centre) for a small number of typical points in M51 (column 4). 
Gas densities a as determined from the CO and H I maps , for 
different values of X: the standard value of 2.8 x 10 20 (column 
5), and 1.1 x 10 20 (column 6). Points Al, A2 and A3 are on the 
arms, II, 12 and 13 are in the interarm regions. 
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